Collective excitation of trapped degenerate Fermi gases 
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We show that the slow driving of a focused laser beam through the cloud of trapped cold fermions 
allows for a creation of the collective excitation in the system. The method, proposed originally by 
us for bosons, seems to be quite feasible experimentally — it requires only a proper change in time 
of the potential in atomic traps, as realized in laboratories already. 
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The experimental realization of the Bose-Einstein con- 
densate (BEC) |l| in a trapped Bose gas has trig- 
gered considerable interest in the field of cold degener- 
ate atomic gases The relatively weak interaction 
between atoms allows for precise experimental manipu- 
lation of an atomic gas. In particular collective excita- 
tions of the BEC such as vortices and solitons (in 
the mean field language - see also ) have been studied 
in details both theoretically and experimentally. Nowa- 
days, it is also possible to cool down Fermi gases to a 
regime where effects of the quantum statistics become 
noticeable [p|-[TTl. Therefore, it would be interesting to 
contrast their collective behaviors p^-p^] with those of 
the BEC. 

Very recently we have proposed a new simple scheme 
that allows for creation of both solitons and vortices in a 
BEC using an appropriate time dependent modification 
(by an additional tightly focused laser beam) of the trap- 
ping potential ||l5|, 161 . The aim of the present work is to 
discuss possible collective excitations of trapped fermions 
and to show that the very same scheme as for bosons 
may be utilized successfully to create excitation in a cold 
Fermi gas. A numerical attempt to create such excita- 
tions using the phase imprinting method JT?! , has been 
recently made [Q. It is not clear, however, which states 
are really excited in this scheme. 

The dominant, at low temperatures, s-wave collisions 
are prohibited for spin-polarized identical fermions due 
to the Pauli exclusion principle. Consequently, it is nat- 
ural to consider first a noninteracting particle model as 
the atoms may interact only through vanishingly small 
p-wave collisions. Our excitation method originally pro- 
posed for noninteracting bosons remains valid even for 
quite strong interaction between particles. We expect 
thus that consideration of the noninteracting Fermi sys- 
tem as a first approximation is perfectly legitimate. 

It is well known that BCS transition may occur in the 
quantum cold gas of fermions at a very low temperature 
Tc ]l^ , p^ . While in the noninteracting particle model 
this effect is absent, we shall show that creation of col- 
lective excitations may occur both for T = (neglecting 
BCS pairing) and for temperatures that are greater than 
Tc but sufficiently low that the effects of quantum statis- 
tics are noticeable. 



The paper is organized as follows. Firstly we specify 
what we consider as a collective excitation in the cold 
fermionic gas on a one-dimensional (ID) example. Sec- 
ondly, we present the basic idea of the excitation scheme 
still for ID noninteracting particle model at zero tem- 
perature limit. Finally we move to the 3D model for an 
experimentally realistic temperature indicating that sig- 
natures of the collective excitations can be observed in a 
laboratory. 

A gas of bosons being in a product state of ground, 
single particle (e.g. mean field) states is a standard ap- 
proximate description of the BEC. It is natural then to 
consider as a collective excitation a situation in which all 
particles are simultaneously transferred to some excited 
state of the mean field effective potential. This results 
in the so called "dark soliton" j20j. Since all particles 
are both initially and finally in the same state, the mean 
field description of the soliton creation requires an ef- 
ficient mechanism for transfer the population from the 
ground to the first excited state. 

The situation is quite different for Fermi system for 
which each state can be at most single occupied (we as- 
sume a spin-polarized identical fermions as realized ex- 
perimentally). Consider a ID Fermi gas at T = pre- 
pared in an ideal collective ground state of the harmonic 
trap. For TV identical fermions this is equivalent to as- 
suming that all oscillator eigenstates, ^„(a;), from n = 
(single particle ground state) to n = — 1 are occupied 
with other levels being necessarily empty. The wave func- 
tion of the N fermion system may be represented then 
by the Slater determinant 
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The corresponding single particle reduced probability 
density reads 
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which (employing Christoffel-Darboux formula |2^]) may 
be reduced to 
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where Hn{x) denote standard Hermite polynomials. 
This density is shown in Fig. ^a. 
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FIG. 1. Single particle reduced probability density for 
A'' — 300 fermions in the ID model, at temperature T — 0, 
before [panel (a) -ground state] and after [panel (b)- collec- 
tively excited state] the sweeping of the perturbation. For x 
we use the harmonic oscillator unit of length, i.e. y/h/muj. 

Now since each fermion occupies a different state, there 
is some ambiguity in defining a collective excitation. We 
define it by requiring that each fermion undergoes a sin- 
gle excitation. Since they are indistinguishable, the final 
wave function is then the Slater determinant involving 
the oscillator functions from n = 1 to n = iV, the corre- 
sponding single particle density is shown in Fig. |l|b. This 
is a most complete analogy to the collective excitation of 
bosons where all N particles gain a single excitation too. 
The difference is that bosons enjoy the same initial (and 
final) state while the fermions need to be excited from 
necessarily different states. Importantly, all bosons pass 
from an even to an odd state in an ideal case, leading to 
the creation of a "dark soliton" , with a node in the cen- 
ter. The collective excitation of a fermionic gas leaves the 
ground state empty that results in a dip in the single par- 
ticle probability density (Fig. |^b). In the ID model the 
relative depth of the dip scales with number of fermions 
as 1/VN. 

How to realize such a final "collective" state? Sur- 
prisingly we show below that a patient way of exciting 
fermions one by one, is quite simple and feasible. More- 
over, it turns out that in a practical implementation it 
seems to be an almost trivial extension of the scheme we 
proposed for an efficient excitation of a Bose gas! 

As in the original suggestion for the bosonic case p^ ] 
we propose to modify the harmonic trapping potential 
by an additional tightly focused laser beam which we 
sweep through the trap. Provided that the frequency of 



this laser beam is sufficiently detuned from an internal 
atomic transition (still being close to the resonance so 
that two atomic levels can be considered only) the upper 
level may be eliminated adiabatically from the analysis 
(see e.g. [^). Then the effective potential for the atomic 
motion, in the ID model, reads 
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where Uo{xo) < (for an appropriately chosen laser de- 
tuning). The amplitude Uo{xo) is proportional to the 
laser intensity, xq = xo{t) is a time-dependent position 
of the center of the laser beam, while a is directly re- 
lated to the cross-section of the (gaussian) laser beam. 
We assume the trapping harmonic oscillator units, i.e. 
ujt for time and yJlifrmL for length, where m stands for 
an atomic mass. 

Assume we prepare the fermionic sample in the state 
(|l|) in a pure harmonic trap. Next we turn on the laser, 
initially focused at the left edge of the trap (a;o(0) ^ 0) 
and move the focus across the trap towards the center 
realizing the potential (^. For all numerical simulations 
presented below we take cr = 0.2 and 



Ufiixa) = (a — fo.To) arctan(a;o). 



(4) 



with a, 6 > 0. In the bosonic case we took \Jq(xq) propor- 
tional to arctan(xo) (ijjl^- In the present case we have 
to add an additional prefactor ensuring that the pertur- 
bation is not negligible in comparison with the harmonic 
part (which is large far from the center of the trap). We 
have chosen the simplest linear form of such a prefactor 
taking a — 18, h — 0.5 in the numerical examples, the 
results are not very sensitive to precise values of a and 
h. During the excitation process xq changes linearly with 
time according to x^it) — x^ifS) -f 0.02i. 

To understand the effect produced by sweeping of the 
laser induced potential across the trap it is sufficient to 
consider the change of single particle energy levels as a 
function of the center of the laser beam xq. The corre- 
sponding level dynamics is depicted in Fig. ^. Observe 
that the levels undergo a series of orderly arranged very 
narrow avoided crossings (not only between low lying 
levels but also between the highest ones). Assume Xq 
is changed sufficiently slowly so as to follow the energy 
levels adiabatically except in the immediate vicinities of 
avoided crossings. Here, since the avoided crossings are 
extremely narrow we assume that they are passed dia- 
batically. Under these premises the excitation scenario 
works as follows. We start with the laser beam situated 
at negative value of xq (with large \xq\) and increase Xq 
up to xo = where the sweeping process ends (when 
the center of the laser beam reaches the center of the 
harmonic trap its intensity is reduced to zero). With in- 
creasing Xq we pass diabatically the consecutive avoided 
crossings encountered on the way. The important avoided 



2 



crossings occur when one of the levels is occupied. As- 
suming that the highest occupied level is A*" — 1, the di- 
abatic passage via the avoided crossing of the N with 
the N — 1 level populates the N level leaving practically 
empty the N — 1 state (the Landau-Zener effect). The 
next crossing, occurring at slightly larger xq value popu- 
lates now the empty — 1 level from N ~2, another one 
populates A^ — 2 at the expense of A^ — 3 and so on. The 
last avoided crossing assumed to be passed diabatically 
populates n = 1 state and leaves empty n = 0. It is clear 
that, provided all A^ avoided crossings are passed diabat- 
ically (all A^ particles are successively excited), we realize 
the desired, described above final Slater determinant in- 
volving the oscillator functions from n = 1 to n = A^ as 
shown in Fig. |l|b. 
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FIG. 2. Energy levels of the potential (g) as a function 
of xq. Panel (a) shows the behavior of low lying energy lev- 
els while in the panel (b) the energy levels around n — 180 
are presented. Note a series of very narrow avoided crossings 
between the neighboring energy levels (for clarity consecutive 
levels are drawn using solid and dashed lines). The inset in 
panel (a) shows a vastly enlarged avoided crossing between 
the ground and the first excited state around xo ~ —6.2. In 
the figure we use the harmonic oscillator units, i.e. hui for 
energy and y^h/mu for length. 



scheme works very efficiently in the bosonic case. There, 
however, a single avoided crossing (between the ground 
and the first excited state) had to be passed diabatically 
by A^ bosons. In the present situation we have to cross 
diabatically A^ different avoided crossings. 

Apart from the experimental approach which is ob- 
viously beyond the scope of the present letter, the other 
way to test the proposed scheme is the numerical integra- 
tion of the time-dependent Schrodinger equation. Since 
we deal with A^ noninteracting particles the problem of 
integration of A^ dimensional Schrodinger equation re- 
duces to the problem of integration of A^ independent 
single particle equations for time evolution of the wave 
functions tpi, see (|^). Such an approach, with the param- 
eters assumed above [i.e. cr, Uo{xo) and xo{t)] yields to 
transfer of particle from n level to n -I- 1 level with the 
probability p higher than 99.95 %, as we have checked for 
n < 180. This indicates that the proposed mechanism of 
excitation is extremely efficient and to the end of this 
paper we will show analytical results assuming p = 1. 

So far we have analyzed the system at the zero tem- 
perature limit. The lowest attainable temperature of 
a trapped Fermi gas in the recent experiments corre- 
sponded to 0.2Tp jl^jl^, where Tp is the Fermi temper- 
ature. To consider the finite temperature effect on the 
excitation process we switch to the 3D model of nonin- 
teracting fermions. We assume a cigar-shaped harmonic 
trap with an additional local gaussian well created by the 
laser beam 
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where we have chosen uj±/u! = 10 as an example. The 
other parameters of the potential are exactly the same as 
in the ID case. The laser beam is directed along z axis. 
Its crosssection in the xy plane is assumed to be tightly 
focused in the x direction having at the same time large 
waist along y. In effect the gaussian character in the 
potential is realized in the x direction only, while in the 
y direction the laser intensity remains constant on the 
size of the Fermi cloud. Changing the trapping potential 
by sweeping of the laser beam allows for excitations in 
the x degree of freedom of the system. The excitation in 
this degree is exactly the same as in the ID case consid- 
ered previously because the 3D model of noninteracting 
particle is separable. 

Assume we have A^ fermions in the harmonic trap dis- 
tributed among the energy levels according to the Fermi- 
Dirac statistics with T > 0. The corresponding single 
particle probability density (integrated over y and z co- 
ordinates) reads 



We have shown llqJlq] that a similar laser sweeping 
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where the sum runs over occupied ipn^ states only. The 
degeneracy factor takes into account that several 
fermions may share the same ipn^ state with different 
other quantum numbers (basically gn^ — N). Next 
we perform the potential sweeping that realizes excita- 
tion in the x degree of freedom, i.e. each occupied ipn^ 
state goes to ipn^+i with a probability practically equal 
to unity. In Fig. ^ we show the results of the excita- 
tion process for N — 40000 fermions in the temperature 
T = 0.3Tp. For example that corresponds to T = 0.2 /iK 
for Lu = 2t: ■ 50 Hz being within the range of current 
experiments ]g|-[Tl|] . The results presented in Fig. ^ were 
obtained by averaging over 100 realizations of N fermions 
in the temperature T — O.STp- A successful excitation 
results in a central dip clearly visible in the figure. Note 
also that in the present case we can employ much more 
particles than in the ID model and the dip in the proba- 
bility density possesses still a considerable depth. That is 
due to the simple fact that before the potential sweeping 
the ipn^ state with — contributes many (go) times to 
the density (^. Consequently, after the excitation, the 
lack of the = state is much more significant than in 
the ID case where such a state could be occupied once 
only. 
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FIG. 3. Single particle reduced probability density af- 
ter the excitation process for A'^ = 40000 fermions. The fig- 
ure represents average over 100 samples corresponding to the 
Fermi-Dirac distribution with T — O.STf- The inset shows 
the same density after a free expansion lasting one period of 
the harmonic trap {2n/ui). In the figure we use the harmonic 
oscillator unit of length, i.e. yjh/muj. 

As a difficulty, however, one might realize that the 
width of the dip is of the order of a few percent of the 
whole distribution of atoms in a harmonic trap. This 
might preclude its direct observation. To overcome this 
difhculty one may, for the detection purposes, supple- 
ment the proposed scheme by the "free fall expansion" . 
Assume that after the sweep by the laser beam and 
the corresponding creation of the collective excitation we 
turn off the trap completely. The group of fermions will 
fall in the gravitational field, the cloud will expand si- 
multaneously. It is important to realize that the cor- 



responding single particle density also expands roughly 
preserving its shape (and in particular the central dip 
structure). The single particle density after the free ex- 
pansion lasting 2tt/uj (i.e. 20 ms for lo = 2tt ■ 50 Hz) 
is shown in the inset of Fig. ^. Observe that the width 
of the dip can quickly reach the size comparable to the 
initial size of the whole cloud. 

To summarize, we have presented a simple method that 
allows for creation of collective excitation in a trapped 
Fermi gas in analogy to a soliton-like state of a Bose- 
Einstcin condensate. The results shown in this work are 
based on the noninteracting particle model. We believe, 
that in analogy to the bosonic case, an interaction be- 
tween fermions (much weaker than in the bosonic case) 
will not affect the efficiency of the proposed scheme. 

The method can be also applied to create collective 
excitation in the so-called Tonks gases (see, e.g., |2^] 
and references therein). Indeed, impenetrable Bosons 
trapped in a quasi ID potential can be modeled by a 
noninteracting Fermi gas. Hence, the excitation scheme 
presented here becomes directly applicable to such a Bose 
system. 
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